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We introduce and formalize the concept of information flux in a many-body register as the influence that the 
dynamics of a specific element receive from any other element of the register. By quantifying the information 
flux in a protocol, we can design the most appropriate initial state of the system and, noticeably, the distribution 
of coupling strengths among the parts of the register itself. The intuitive nature of this tool and its flexibility, 
which allow for easily manageable numerical approaches when analytic expressions are not straightforward, 
are greatly useful in interacting many-body systems such as quantum spin chains. We illustrate the use of this 
concept in quantum cloning and quantum state transfer and we also sketch its extension to non-unitary dynamics. 

PACS numbers: 03.67.Hk, 75.10.Pq, 05.50.+q 



I. INTRODUCTION 

Very recently, we have witnessed a strong interest of the 
quantum information processing (QIP) community in prob- 
lems of interacting many -body systems [1]. Such interest is 
motivated by the unresolved necessity of tackling, theoreti- 
cally and experimentally, quantum processors equipped with 
large registers. The investigation of multipartite devices has 
been boosted by the observation that specific forms of built- 
in and permanent intra-register couplings can be used for the 
purposes of quantum communication and quantum computa- 
tion yj, |2[] . Indeed, it has been shown that the control over 
such systems can be sensibly reduced in a way to avoid the 
generally demanding fast and accurate inter-qubit switching 
and gating (we will name such a scenario as control-limited). 
However, the price to pay for the performance of efficient 
operations is the pre-engineering of appropriate patterns of 
couplings. Determining the exact distribution of coupling 
strengths for a given inter-qubit interaction model in a control- 
limited setting is often a matter of craftsmanship or the result 
of the exploitation of peculiar topological properties of the 
physical system at hand yj]. A more systematic approach to 
such a problem is evidently in order. 

Here, we provide a tool for this task by introducing the con- 
cept of information flux in a quantum mechanical system. An 
intuitive understanding of our idea can be given in terms of the 
influences that the dynamics of a selected element of a multi- 
partite register experience due to interaction channels with the 
other parties. We show that a desired protocol is performed by 
arranging the network of interactions in a way to privilege or 
repress specific interaction channels. The realization of an op- 
timal QIP task is therefore translated into the maximization of 
the information flux associated with such channels. This re- 
sult is twofold interesting: conceptually, the study of informa- 
tion flux in multipartite schemes for QIP allows for a deeper 
insight into the physical mechanisms behind an investigated 
protocol. We provide a significant example, in this paper, by 
addressing the case of universal quantum cloning machines 
where the achievement of the optimal cloning fidelity is shown 
to correspond to an information-flux maximization problem. 
Pragmatically, as anticipated, it helps in designing the most 



appropriate pattern of interaction strengths and/or initial state 
of a register. We address the case of information transfer in 
spin chains as a scenario where such a possibility is partic- 
ularly useful. Although our attention is focused on registers 
of qubits, it is important to notice that the dimensionality of 
the Hilbert space of the elements of a register does not rep- 
resent a limitation to the effectiveness of our approach. We 
also sketch the way interaction channels can be generalized to 
the non-unitary case by explicitly introducing noise and de- 
coherence in the system and considering the corresponding 
open-dynamics information flux. 

This paper is organized as follows. In Sec.|II]we introduce 
the approach we use to analyze different schemes for QIP. 
In Sec. [HI] we investigate the quantum circuit proposed by 
Buzek et al. Q4£] for an optimal two-copy universal quantum 
doner (UQCM 1 — * 2). An implementation of this protocol, 
by means of spin chains, is studied in Sec. |IV] where a clear 
picture of the advantages offered by the information-flux ap- 
proach is provided. We consider quantum state transfer across 
a spin chain in Sec.[V] where we also show how to adapt our 
approach to handy numerical analysis which is used to reveal 
interesting features of long and disordered chains. In Sec.lVIl 
we extend the information-flux concept to non-unitary dynam- 
ics. Finally, in Sec. lVIII we summarize our results. 



II. CONCEPT 

Let us consider a register of N interacting qubits coupled 
via the Hamiltonian 7is g \{t") whose structure we do not need 
to specify here. Our assumption is that 7i{ ff j(t) depends on 
a set of parameters gj (which could stand for the coupling 
strengths between the elements of the register) and a gener- 
alized time parameter t. We adopt the notation according to 
which Ej = ®iz\i k <g> <7 Ej ®iL j+l h (S = X,Y,Z,I) is 
the operator that applies the er^ Pauli-matrix only to the j-th 
qubit of the register. Here, aj. = ij with Ij the 2x2 iden- 
tity matrix of qubit j. In the remainder of this paper, unless 
explicitly specified, we work in the Heisenberg picture where 

time-evolved operators are indicated as = WY>jU with 
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FIG. 1: (Color online) Scheme of the processes we consider. A com- 
putation or communication step is interpreted as a black box (whose 
operation depends on the coupling scheme within a multipartite reg- 
ister of qubits) with movable input and a detection terminal. Through 
the study of information-flux dynamics, we can design the best cou- 
pling scheme for a chosen QIP operation. 



U(t) = exp[— (i/h) J Ti{ g }(t')dt']. We say that there is in- 
formation extractable from qubit j at time t whenever there 
is at least one £ (obviously excluding the identity) for which 

(T,j(t)) 0. Here, the expectation value is calculated over 
the initial state of the register |\Prj)i m- In some cases, the 
analysis of expectation values in the Heisenberg picture can 
be considerably helpful in gathering information about the dy- 
namics of a system [5[]. 

We adopt, hereafter, the following schematic description of 
a computation or communication process: we suppose to have 
access to a selected qubit of a multipartite register and we con- 
sider it as the input terminal of the black box given by the 
rest of the elements. We then use a detection stage which can 
be attached to a suitable output port, connected to one of the 
qubits in the black box (a sketch is given in Fig. [TJ). In this 
picture, the initial state of a quantum system is described by 
the state vector \^o)i m = \4>o)i IV'o^ jy This * s ^ e case 
in which the first qubit is initialized in a generic input state 
(and is separable with respect to the rest of the register) and 
^0)2 n represents the initial state of the other qubits that, in 
general, can be mutually entangled. We assume this state to 
be known and independent of the input state. The assumption 
of a known initial state of the register {2, .., N} is physically 
motivated as it corresponds to the situation assumed in many 
control-limited QIP protocols HUSH]. We can thus inter- 
pret a quantum process as the flux of appropriately processed 
information from the input qubit to the remaining components 
of the register. Such a flux is witnessed by any explicit depen- 
dence of the dynamics of the i-th qubit from the operators 
associated with the input one. Therefore, in order to find if 
qubit i has developed any extractable information at time t, as 
a result of an information flux from the input qubit, we need 

to study the dependence of {^Po 1 (i) | ^0) ' s on at least one of 

(\Po|5^i(0) I^Po)^. In what follows, we show that the study 
of the dynamics of these expectation values, taking advan- 
tage from the possibility of working with "known and fixed" 
states, allows for proper interaction engineering and optimiza- 
tion that can be used for variegated tasks of control-limited 
QIP. 

Each T,i(t) can be decomposed over the operator-basis 



built out of all the possible tensorial products of single-qubit 
operators acting on the elements of the system {1,..,N}. 
Such a basis has dimension 4^ and its elements can be par- 
titioned in four disjoint groups, each of dimension 4 Ar ~ 1 . 
Each group contains operators having the form <8> Gk,2..N 
(k = 1, ..,4 Ar ~ 1 ) with Gk,2..N being the tensorial product of 
single-qubit operators acting on the elements of the system 
{2, .., TV}. By using the knowledge we have about the initial 
state of the {2, .., N} register, we introduce the vector x suc h 
thatxfe = 2..iv(^o|^fe,2..JV |"0o) 2 ..jv Therefore 



(#o|Si(*)l*o) 



=X,Y.Z,I fc=l 



(1) 

Here, the matrix a ES (t) incorporates the details of the time- 
evolved multi-qubit operators and is determined once Ti{ g } (t) 

is assigned. The coefficient Xf s ' (t) = J2t=i a ife E 'MXfc 
defines and quantifies the flux of extractable information (or, 
shortly, the information flux) from Sj to Sj at time t. The 
formal definition of the information flux highlights the dual 
nature of the control we can operate over the dynamics of 
a multi-qubit system. Indeed, besides the dynamical part of 
2p E (t), there is a time-independent part that incorporates in- 
formation about the initial state of the register, input qubit 
apart. This represents an additional control over the dynam- 
ics at hand. Information flux towards a specific element of 
the register can be depleted or enhanced by properly prepar- 
ing the register {2, .., N} in a way to engineer Xfc' s - Notice 
that the intrinsic dependence of the Hamiltonian on the set 
{<?} makes the information flux implicitly dependent on the 
coupling scheme being chosen. Obviously, by choosing a dif- 
ferent partition for the four groups of operators in which the 
basis for an iV-qubit system has been divided, Eq. ([TJ) can be 
straightforwardly modified in order to make explicit which is 
the information flux from Sj to ££. In our black-box model, 
this corresponds to a change in the position of the input ter- 
minal (cf. Fig. [TJ. The dual control discussed above can be 
fully utilized in the preparation of a multipartite device in the 
most appropriate configuration (of couplings and initial state) 
for a given QIP task. The potentialities of this approach are 
better illustrated by means of explicit examples. In what fol- 
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FIG. 2: (Color online) (a): Logical circuit for optimal UQCM 
1 —* 2 1 4]. Each line represents a qubit and time flows from left 
to right. We show single-qubit rotations, indicated by the letter R 
inside a circle, and two-qubit CNOT gates, (b): Information flux 
from the input qubit to the rest of the register in the copying stage of 
the UQCM 1 — f 2. The yellow arrows represent the a z information 
fluxes and the blue arrows represent the a x information fluxes. 
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lows, we address the important problems of quantum cloning 
and quantum information transfer. The use of the informa- 
tion flux approach in the investigation of the latter problem, 
in particular, turns out to be quite useful for the analysis and 
characterization of quantum state transfer in a long spin chain 
with imperfections. 



in. UNIVERSAL QUANTUM CLONING MACHINE 1 -> 2 

We consider the problem of an optimal two-copy univer- 
sal quantum copy machine (UQCM 1 — + 2), for which ex- 
tensive literature is available (TO. In our study we concen- 
trate on the quantum circuit proposed by Buzek et al. [4] 
and shown in Fig. [2] (a). In order to adapt the framework 
of Ref. [4] to the notation used so far, we say that qubits 
2 and 3 embody the register described in the previous sec- 
tion and qubit 1 is the input qubit. The circuit is composed 
of two different stages: the preparation stage, composed of 
single-qubit rotations i?'s and controlled-NOT (CNOT) gates, 
prepares the second and the third qubit in the state ^0)23 = 

(1/V6) (2 1 00)23 + 1 01) 23 + 1 10)23 ) while the state to be c l oned 
is encoded in the input qubit. The copying stage in Fig. [2] (a) 
produces two optimal clones of the input state on the second 
and the third qubit. The cloning fidelity is |, which is the 
maximum for a UQCM 1 — > 2 jSQ. 

Let us now use an information-flux approach to the prob- 
lem by studying the flux from £1 to E 2 3 . Of course we need 
to investigate only the copying stage, because no interaction 
involves the first qubit during the preparation stage, which can 
be performed off-line. The copying stage is built out of four 
CNOT gates, which altogether give the time-evolution oper- 
ator. We only need to consider X.2,3 (t) and Z 2j3 (£), for in- 
stance, as Y~2^(t) is found as iX 2 , 3(^)^2, z{t)- For a general 
CNOT gate acting on its control and target qubits H, it is 
immediate to get the following relations: 

a c = X c X t , X t — Xti Z c — Z c , Z t = Z c Z t , (2) 

with X c , Z c (X t , Z t ) the operators of the control (target) qubit 

before the action of the gate and X c , Z c (X t , Z t ) the anal- 
ogous operators after the action of the gate. Obviously, in 
this example, the generalized "time-parameter" corresponds 
to the "logical application of a single CNOT gate." By using 
these results, we straightforwardly find the expressions of the 
evolved Xi and Zi (i = 1, 2, 3) for the register of the UQCM 
1^2 after each step of the copying stage. These are summa- 
rized in Table U [ 11]. Our first observation is that there is no 
information flux from Xi to X 2 ^ during the first interaction, 

as X 2 3(ti) do not depend on X\. The same holds for Zj,(ti). 
On the other hand, by applying the formalism introduced 
in the previous section, an easy calculation let us see that 
<*o|^ 2 (*i)|*o) = (MiMZiZzWaWo) = KMZM, 
i.e., X 2 Z (t\) — 2/3. After t\, no information flows from 
Zx to Z 2 , as it can be seen from Table [I] so that if z (i 4 ) = 
2jf (ti). An analogous reasoning leads us to conclude that 



given by Xgf (i 4 



i$J(ti) = ZP(i 4 ) = 2/3. In this 
case, where all Xff 's are equal and Xff 's vanish for £ 7^ 
we can affirm that the Bloch vectors of the output states are 
proportional to the one associated with the input state with the 
shrinking factorizing equal to the information flux. The tim- 
ing and pattern of information fluxes are pictorially illustrated 
in Fig. |2](b). The specific pattern of fluxes is a particular re- 
sult of the symmetry of the UQCM 1 — > 2. In general, when 
the flux is restricted to homonymous operators, the following 
simple relation can be derived connecting the fidelity Tu (t) 
between the input state and the state of qubit i at time t and 
the corresponding information fluxes in the register 



T,=X,Y,Z 



(3) 



In the specific case at hand and regardless of the input state, 
•^"12(^4) = Tizitt) = 5/6, as it should be. Equation (O 
can be straightforwardly generalized to the case where cross- 
operator fluxes are present. Such an expression is of no rele- 
vance, however, for the purposes of this paper. 

The pedagogical example discussed above, however, does 
not exhaust the possibilities offered by an information-flux 
approach. Indeed, as anticipated, the method can also be 
used to find the optimal preparation stage of a circuit de- 
signed to perform a given task. If we consider the above- 
mentioned copying stage, without knowledge of the input 
state, we can use the information-flux analysis to design the 
proper preparation for an optimal UQCM 1^2. Let us 
suppose that, starting only from the information-flux anal- 
ysis previously done on the copying stage, we want to de- 
duce the initial state of the second and the third qubit to ob- 
tain an optimal UQCM 1^2. The conditions on such an 
initial state can be easily obtained by imposing that the two 
clones are identical (symmetry condition) and the cloning pro- 
cess is independent of the input state (universality condition). 
The symmetry condition implies that Xf^(t4) = ^J" 2 ^) 
while the universality condition gives 2fT (t 4 ) = Xfrp (t 4 ) 
(with E,E' = X, Y, Z). From the study of the copying 
stage we know that no cross-operator flux is possible itllll . 
so that the use of Eq. (f3]l is legitimate. For the generic input 
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the remaining non-zero information fluxes in this problem are 



TABLE I: The evolved operators after time tj [i.e., after the applica- 
tion of the j-th CNOT in Fig.f2](a)] expressed in terms of the opera- 
tors before the copying stage. 
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«0o| & Xl I0o}) 2 + «0o| fiy, |0o)) 2 + «0o| & ?1 |0o)) 2 = 1. 
As the process must be optimal, the fidelities have to be 
maximized (with the constraints imposed on the information 
fluxes), regardless of the state to be cloned. This means the 
maximization of the information flux. For a generic input 
state |Vo> 23 = a|00) 23 + /3|01) 23 + T |10) 23 + <J|H) 23 , 
this requirement implies the maximization of (t/'ol-^lV'o) = 
(0o|Z 2 |-0o) = (0o|^3|^o) = (-00 1 -£3 1 ^o) which results in 
a = y/2/3, fi = 7 = 1/ \/6, and 6 = 0. We therefore re- 
trieve, through an information-flux analysis, the initial state 
of the optimal UQCM 1 -► 2. 

As an example of how the initial preparation of the reg- 
ister can help in designing the optimal operation to perform 
by implicitly "inhibiting" information flux, let us consider the 
following simple but yet illustrative case. Suppose that, us- 
ing again the copying stage of Figs. [2] we want to obtain 
a circuit that generates the output state of qubit 2 equal to 
the input state, regardless of the third qubit (in some sort 
of fully-biased asymmetric cloning). Again, we want to de- 
sign an appropriate preparation stage for qubits 2 and 3. The 
complete asymmetry of the process requires (quite intuitively) 
that I^{t^) = 1, By using Table U again, this turns out to 
be equivalent to (ipa\^3\ipa) — (V'ol^lV'o) = 1. For the 
generic input state ^0)23 use( ^ above, these conditions are 
satisfied only for a = /3 = 1/ v2 and 7 = 5 = 0, so that 
1^0)23 — (l/v 7 ^) |0) 2 (|0) 3 + |1) 3 ). The preparation stage 
should thus be a simple Hadamard transform over the third 
qubit. 

IV. UQCM 1 -> 2 WITH SPIN CHAINS 

The information-flux approach offers remarkable advan- 
tages also in the cases where, rather than abstract quantum 
circuits, we have to deal with physical interaction models 
described by coupling Hamiltonians. This situation is par- 
ticularly relevant for scenarios of control-limited QIP, where 
multi-qubit interactions are often employed for the engineer- 
ing of appropriate quantum interference proc esses resulting in 
the performance of QIP tasks |Sl HH 111 Q. In particu- 
lar, it has been realized that information transfer and quantum 
cloning can be performed through well-known multi-qubit 
models such as XY or Heisenberg couplings |lj, |fj. Here, 
we focus our attention to the optimal UQCM 1 — * 2 real- 
ized with a three-qubit chain [10] whose coupling Hamilto- 
nian, within a nearest-neighbor anisotropic Heisenberg model, 
reads H = {J /2)Y J 2 l=1 {X l X l+1 + Y l Y l+1 +\Z l Z l+1 ). Here, 
J is the coupling strength and A > is the anisotropy parame- 
ter. The central qubit (labeled by 2) is the input one, encoding 
the state to be cloned, while the external qubits 1 and 3 are the 
output ones. Note that this is a situation requiring the input 
terminal in our black-box schematics to be shifted to qubit 2. 
In an information-flux perspective we will thus be interested 
in the fluxes from qubit 2, this time. Here, we demonstrate, 
by means of the proposed new tools, that for A = 2, a proper 
initialization of qubits 1 and 3 and a careful timing allow for 
optimal UQCM 1 -> 2 fioll. 

The symmetry condition, in this case, is a direct conse- 



quence of the physical system at hand so that, in evaluat- 
ing the information flux from qubit 2 during the evolution, 
we can just concentrate on qubit 1. From the analysis con- 
ducted so far on UQCM 1 — *• 2, it appears that the conditions 
of no cross-flux and 2f;f(t*) = 2/3 (V S = X, Y, Z and 
at a proper time t*) guarantee the universality condition and 
the optimal value of fidelity to be achieved. A detailed cal- 
culation shows that the information fluxes between homony- 
mous operators of qubits 2 and 1 (or equivalently 3) depend 

on 13 (V'o I 100)13 w ^ tn IV*) 13 a generic two-qubit pure 

state. In order for such expectation values to be physically 
meaningful when we impose Tff(t*) = 2/3 (i.e. in order to 

have 13 (00 1 ^1^3 \ipo} 13 € [— 1, 1])> the only possible choice 
for the anisotropy parameter is A = 2 11311 . In turn, it im- 
plies that 13 (V>o I ^1^3 IV'ohs = 13 (0o I Y x % ^0)13 = 1 and 
13 (0o I Z1Z3 |^o)i 3 = — 1- The only state of qubits 1 and 
3 satisfying these conditions is |0+) 13 = (l/v2)(|01) + 
1 10) ) 13, so that |0>o)i3 = l0+)i3- With these choices, 
Zf* (<) = ( 2 / 3 ) sin 2 (V3Jt) and at t* = tt/(2V3J) we have 
the optimal fidelity T 2 i{t*) = T 2 ?,{t*) = 5/6. Our rigorous 
result confirms the "guess" for the initial state of the clones in 
Ref. II 1011 and demonstrates the uniqueness of such a prepara- 
tion for the optimality of the process. 

There is a profound difference between the cloning per- 
formed through the multi-spin coupling and the formal cir- 
cuital dynamics depicted in Fig.[2](a), which the information- 
flux approach is able to clearly spot out. For cloning realized 
with the Heisenberg coupling, the information fluxes from the 
input to the output qubits are always equal during the evolu- 
tion. As shown in Fig. [2] (b), on the other hand, the quan- 
tum circuit allows the information flow in the register in a se- 
quential manner. Therefore, while the universality condition 
is always respected during the UQCM 1 — * 2 realized with a 
multi-spin coupling, this condition holds in the quantum cir- 
cuit only considering the whole transformation, but not during 
the single steps performed by the quantum gates. The same 
holds for the symmetry condition. Therefore, if we stop the 
interaction in the spin system at a time different from t*, the 
cloning fidelity will certainly be lower than the optimum but 
it will not depend on the input state and it will be the same for 
both the clones. Our method can also be applied to the very 
recently proposed case of UQCM 1 — > A" through spin-chain 
models llHl . 



V. STATE TRANSFER 

An important part of our study is the investigation of quan- 
tum state transfer (QST) from the information-flux viewpoint. 
We show how such a change of p ersp ective with respect to the 
standard approach to QST lllll3L ll2lll5ll . allows for the study 
of long-chain dynamics in the presence of imperfections in a 
convenient way. Let us gradually arrive to such a result. 

It is well-known that perfect QST is possible for a three- 
qubit chain of equal XY couplings governed by H = 
(J/2) Y^ =1 (XiX l+1 + YiY l+1 ) yD. This can be seen by ex- 
plicit calculation of the spectrum of Tt and the corresponding 



fidelity. In our analysis, the same result is achieved simply 
by noticing that, in these conditions, l xx (t) = T YY (t) — 
\/T^ z (t) = sin 2 (J</\/2), so that at t* = w/V2J there is a 
unit information-flux from each operator of the first qubit to 
the homonymous operators of the last one. Incidentally, we 
notice that an imprecision in t* affects X 3 z (t) much more 
than the other fluxes. This means that the transfer fidelity 
in a process which is stopped at the wrong interaction time 
drops faster for input states with a large modulus of the a z - 
expectation value. This explains the frailness of the QST pro- 
cess with input states having a large population-inversion. 

The achievement of perfect QST by using a long linear 
chain of qubits all prepared in their ground state is known to 
be possible for weighted coupling strengths following the pat- 
tern Ji = Xy / i(N — i) with i = 1, .., N — 1 labeling a site in 
the chain: After a time t* = j , we obtain perfect state transfer 
from the first to the last qubit [3]. While this result can be re- 
trieved through the information-flux approach, here we would 
like to be more pragmatic. The achievement of a unit transfer 
fidelity, although obviously desirable, might not be necessar- 
ily in order. Depending on the specific protocol that has to be 
realized after the QST process, a lower threshold to Tin may 
be enough (if Tin is sufficiently large, for instance, one can 
use state-purification procedures after the receipt of the state). 
This can be helpful in practical QST implementations if, by 
paying the price of a lower transfer fidelity, the demanding 
requirements over the coupling pattern are relaxed. Here we 
show that this possibility is actually realistic. 

For some specific coupling configurations and assuming pe- 
riodic boundary conditions, the XY Hamiltonian can be diag- 
onalized by means of a sequence comprising Wigner- Jordan, 
Fourier and Bogoliubov transformations [16j]. However, this 
does not hold for any pattern of coupling rates. As here we 
consider an inhomogeneous set of J^'s in an open chain ruled 

by the Hamiltonian H = \Y,fji HWi+i + 
we decide not to rely on these methods [17]]. By means of a 
simple application of the operator-expansion formula, on the 
other hand, it is possible to design an efficient numerical appa- 
ratus to analyze the information flux. The method is based on 
the construction of recurrence formulas building up the flux 
of information between operators of specific qubits along the 
chain. For instance, for the non-trivial case of N — 5 we have 

X 5 (t) = a(t)X 5 + b(t)Z 5 Y 4 + c(t)Z 5 Z 4 X 3 
+ d(t)Z 5 Z 4 Z 3 Y 2 + e(i)^iWi. 
The time-dependent coefficient ait) can be obtained by the 
formula a(t) = $Z" a 2i(— 1)' tjoti , where a* is related to 
bi, Ci, di, and through 

<H = Jk&i-l, h = J4CI1-1 + J3C1-1, 
Ci = hh-i + Jidi-i, di = J 2 Cj_i + Jiej_i, (5) 
= J\di-\ 

with the initial conditions <xq = 1, 60 = Co = do = cq = 
0. The other coefficients b(t), c(t), d(t) and e(t) are sim- 
ilarly defined. As the initial state of the chain is \^f ) = 
|0o)i |0000> 2 5 , we have x (t) = e{t). 
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(a) (b) 




FIG. 3: (Color online) (a): l£ x in a spin chain of length N = 101 
with equal central coupling strengths, against the rescaled interaction 
time Jt and 77. (b): T xx for the same conditions in panel (a) and 
Tj — rj m ax, plotted against Jt and the standard deviation a of the 
Gaussian distribution which determines the disorder parameters 5i. 



Having illustrated the formal technique, we now move to 
the case of an arbitrary number of qubits. From now on, our 
limited-control assumption is that we are able to engineer the 
strength of the coupling rates of just the extremal qubits ( J\ 
and J N -i) [18]. We take J t = J (for i = 2,.,N - 2), 
Ji = Jn-i = T}J and study the behavior of the flux 
against the dimensionless interaction time Jt (within a rea- 
sonably long range) and the inhomogeneity parameter 77. It 
is important to stress that the limitations to our analysis are 
simply the results of the work conditions we have assumed. 
In principle, any pattern of coupling strengths can be intro- 
duced and studied through our approach. In Fig. [3] (a) we 
show 1 XX for the case of N = 101. Evidently, there are val- 
ues of (77, Jt) for which the information flux towards the last 
qubit in the chain is almost ideal. For instance, at 77 ~ 0.5 and 
Jt ~ 27.6, Iiq X — 0.93. An explicit evaluation reveals that 
lJj Y — Tn X while %n Z \ being equal to their product, reaches 
its maximum value for the same 77 and rescaled time as X x x . 
This implies that a transfer fidelity > 0.865 is achieved for 
the worst-case scenario given by the fully polarized state |1) 1 
being transferred across the chain (here we have used the gen- 
eralized formula for fidelity, due to the fact that cross-operator 
fluxes are present). To give a complete picture of the way the 
transfer dynamics behaves in such a control-limited scenario, 
in Figs. [4] (a), (h) and (c) we present the maximum value of 
I xx and the corresponding Jt and 7/ against the length of 
the chain up to N — 101. Compared to the perfect transfer 
case of N = 3, which occurs for all-equal coupling rates, a 
decrease down to J\ = Jn-i = J/2 is sufficient to achieve 
high-fidelity QST. 



(a) (b) (c) 




FIG. 4: [(a), (b), (c)] Maximum value of the information flux T x , 
the corresponding rescaled time Jt and inhomogeneity parameter 77 
against the number of spins N for chains of up to iV = 101 qubits. 
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The flexibility of our approach is well-witnessed by the eas- 
iness with which imperfections in the coupling rates' pattern 
can be incorporated. This is an important point in any real- 
istic assessment of QST as it will be difficult, in practice, to 
arrange an ideal distribution of coupling strengths. We have 
therefore considered disordered Ji's, each deviating from the 
value corresponding to the QST described above by a quantity 
8i that follows a Gaussian distribution centered at and with 
a standard deviation a of, at most, 5% Ji lfl9ll . For N = 101 
qubits, the corresponding is shown in Fig. 0(b). The 
performance of the chain is quite robust against the disorder. 
In particular, an interesting feature that can be experimentally 
helpful is that the time at which the flux is maximized does 
not change significantly with the strength of the disorder with 
the maximum being only slightly dependent on it. 

VL OPEN-SYSTEM DYNAMICS 

Here, we sketch how the framework presented in this pa- 
per can be extended so as to address the case of open-system 
dynamics where the register is subject to inter-element inter- 
actions and is also exposed to environmental effects. In this 
case, as we explain, the information flux will contain a purely 
coherent part, arising from the unitary dynamics of the system 
at hand, and an incoherent contribution. The starting point of 
our study is a master equation (ME) for the open dynamics of 
the register. For the sake of simplicity, we assume a Marko- 
vian environment, described by a Liouvillian C(g) acting on 
g, the density matrix of the system. To fix the ideas, we con- 
sider the case where each qubit interacts with its own environ- 
ment. The generalization of this situation is discussed later in 
this section. We thus have 

i L 1 (6) 

~|[o-^,[d- z ,i, (?]]}• 

Here, <7+,; (<X-,j) is the raising (lowering) operator of qubit 
i, r is the rate of dissipation of each qubits immersed in a 
zero-temperature dissipative environment (i.e., a thermal bath 
of mean occupation number n = 0). Finally 7 is a dephasing 
rate. We assume that the register undergoes a simple free evo- 
lution and we move to a proper interaction picture, so that the 
register's ME is d t g = C(g). We multiply this by the time- 
independent operator of qubit k and trace over the register 
so that the ME can be transformed into the set of Langevin- 
type equations (one for each £&) 

d t (h(t)) = J2 ^~(0+At)h(t)a-,i(t)) 

- (E k (t)a +t i(t)a-,i(t)) - (o- +li (t)o-_,i(t)E fc (t))) 
-7«£fc(*)> - &,i(*)£k(*)£*,i(*)>)} > 

(7) 

where the expectation values are all calculated over the initial 
state of the register. Obviously, this expression can be easily 



generalized to the case of n ^ by introducing an additional 
term to the Liouvillian and to the case where an arbitrary inter- 
action Hamiltonian 7i{ g }(i) is considered. The time-evolved 
operators &± t i(t) and a Zi i(t) can be expressed, in general, in 
terms of sums of single-qubit operators, weighted by proper 
time-dependent functions. We therefore immediately recog- 
nize that Eq. (0 (or its generalization) provides the open- 
system dynamical equations for the information fluxes within 
the register. In particular, the case at hand addresses just the 
incoherent evolution of the information fluxes (equivalent to 
the noise-part of Langevin-type equations for the expectation 
values of the operators of the system). The described ap- 
proach can be applied, in exactly the same manner, to the 
case where the individual-environment assumption does not 
hold. By studying the information flux between particular 
states of a register of qubits collectively interacting with an 
environment, we gather information about the rate of creation 
of environmentally -mediated entanglement and the effect that 
leakage has on it, for instance [20]. Non-Markovian ME, ad- 
dressing the case of environment with memory, can also be 
treated with our approach. 

The use of expectation values over properly designed initial 
states, once more, can be considered a useful simplification in 
the analysis of the dynamics of a register. The application of 
this framework to pragmatic examples, with particular atten- 
tion to collective or memory-preserving environments, goes 
beyond the scope of this work and will be the subject of fur- 
ther investigations. 

VII. REMARKS 

The possibilities offered by control-limited QIP with multi- 
qubit interactions require the development of exploitable ap- 
paratuses for the design of the coupling distributions and 
preparation stages appropriate for a set task. This is manda- 
tory in order to achieve the optimal pattern of quantum inter- 
ferences at the basis of any quantum communication and com- 
putation protocol. Here, we have introduced the information 
flux as a promising tool in this respect. We have demonstrated 
how it can be used in order to gather important insight in the 
performance of existing protocols for QIP and for their devel- 
opment towards optimality. Particular attention has been paid 
to the important problem of optimal cloning in interacting spin 
systems. On the other hand, its intrinsically operative nature 
allows for the design of quite manageable numerical recipes, 
particularly important in the manipulation of very large com- 
putational registers. These have been instrumental in the study 
of a long chain of coupled qubits acting as a support for quan- 
tum state transfer. We have shown that a simple coupling con- 
figuration can be designed allowing for long-haul communi- 
cation, resilient to (quite large) local fluctuations in the cou- 
pling strengths. Finally, we have sketched the way the concept 
of information flux can be extended to open-system dynamics 
by means of Langevin-type equations of motion. We believe 
our technique has rather promising possibilities of application 
in problems of many-body physics in the presence of disorder 
and decoherence lull . 
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